Gamma rays from star-forming regions 
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Abstract. Star-forming regions have been tentatively associated with gamma-ray sources since the early days of the COS 
B satellite. After the Compton Gamma-Ray Observatory, the statistical evidence for such an association has became 
overwhelming. Recent results from Cherenkov telescopes indicate that some high-energy sources are produced in regions 
of active star formation like Cygnus OB2 and Westerlund 2. In this paper I will briefly review what kind of stellar objects can 
produce gamma-ray emission in star-forming regions and I will suggest that the formation process of massive stars could in 
principle result in the production of observable gamma rays. 
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INTRODUCTION 

The identification of star-forming regions (SFRs) with 
gamma-ray sources has a rather long story that can be 
traced till the first attempts to make statistical correla- 
tion analysis between COS-B data on discrete gamma- 
ray sources and SFRs (e.g. [1|). With the advent of 
the Compton Gamma-Ray Observatory and the Ener- 
getic Gamma-Ray Experiment Telescope (EGRET) the 
statistical evidence for a physical association became 
overwhelming (|2j |3]). In the early 1980s, theories of 
cosmic-ray acceleration in the strong winds of early- 
type stars were worked out in great detail and massive 
stars were proposed as potential gamma-ray sources (e.g. 
|ffl[5]|6][71)- However, other types of sources could con- 
tribute to the gamma-ray emission detected in the di- 
rection of SFRs, sources such as accreting black holes, 
pulsars, supernova remnants, and pulsar wind nebulae. It 
was not until the discovery of the extended unidentified 
source TeV J2032+4130 in Cygnus OB2 by the HEGRA 
array of imaging Cherenkov telescopes in 2002 (8) that 
the hypothesis of a high-energy source powered by stel- 
lar effects in a SFR was considered as sound from an 
observational point of view. Since then, other two high- 
energy sources have been detected in SFRs: Westerlund 2 
[9 1 and W43 HO). All this points out in the direction that 
early-type stars and the collective effects of their winds 
can accelerate particles beyond TeV energies and pro- 
duce significant gamma-ray emission. In this review, I 
will discuss how gamma rays can be produced by stars 
in the different evolutionary stages of a SFR. But before, 
some account on the current ideas of star formation is in 
order. 



THE FORMATION OF MASSIVE STARS 

The formation and early evolution of massive stars are 
not yet well understood. Two different mechanisms have 
been proposed. These mechanisms do not necessarily 
exclude each other. They are the accretion IfTTl [121 and 
the coalescence 11311 massive star formation processes. 
In the first case the massive star is formed by accretion 
of gas in a dense molecular core. Such an accretion 
episode is always accompanied by ejection of thermal 
jets that can propagate inside the molecular cloud or even 
break out from it (for a review see, e.g., IfMlD . In the 
coalescence picture the high-mass stars are formed by 
merging of low and intermediate mass protostars in the 
dense cluster environment. In either case the formation 
takes place inside a dense and giant molecular cloud. The 
SFR, then, starts as a giant dark cloud that is soon heated 
and illuminated by the star formation that goes on in its 
interior. As time goes by, the massive stars are formed 
and their strong winds sweep the ambient material giving 
rise to a mature open cluster. Both stages (initial and 
final) are illustrated in Fig.[T] 

The dense molecular clouds that are the cradles of 
massive stars have relative small sizes (< 0.1 pc), high 
temperatures that make of them strong IR sources (T > 
100 K), densities of the order of ~ 10 7 cm~ 3 or more, 
and high luminosities in the range of 10 4 - 10 6 L .The 
conditions of the medium are surely not homogeneous, 
with a power-aw density distribution n °< r~ p , with p in 
the range between 1.3 and 1.8 03] ■ The massive stars are 
formed under very high infall mass rates of the order of 
10~ 3 — 10~ 2 M yr _1 . The whole formation process is 
rather rapid, with timescales of the order of a few kyr. 

One way to test the different scenarios proposed for 




FIGURE 1. Two different stages of a SFR evolution: a dark 
cloud and a fully formed open cluster. 

massive star formation is the detection of outflows re- 
lated to the accretion process in massive protostars. As 
we will see, there are some prominent examples of such 
outflows. 



COLLECTIVE EFFECTS OF STELLAR 
WINDS IN STAR FORMING REGIONS 

Once a massive stellar association is formed, the com- 
bined effects of the stellar winds are expected to re- 
sult in particle acceleration up to relativistic energies 
|[T6l [T71 [181 . In addition, supernova explosions of very 
massive (and hence quickly evolving) stars can help to 
create super-bubbles and collective shocks where parti- 
cles can be efficiently accelerated (e.g. [19|). Hence, the 
average cosmic ray (CR) density in SFRs is expected to 
be higher than the galactic average. The locally generated 
CRs can interact with both extended passive targets, like 
molecular clouds or the rich intracluster medium (e.g. 
E0II2TII22I ). or with the very same winds of the massive 
stars ll23ll . In the latter case, a modulation effect due to 
convection is expected, and hence the low-energy emis- 
sion would be greatly diminished or suppressed. Such a 
picture might correspond to that of extended TeV sources 
without significant correlated emission at low energies 
like TeV J2032+4130 B4ll25ll26l. The same mechanism 
might operate in starburst galaxies where combined ef- 
fect of the large number massive stars should produce a 
high density of CRs in a very rich environment Il27l l28l . 



GAMMA RAYS FROM MASSIVE STARS 
AND MASSIVE BINARY SYSTEMS 

Individual hot stars with strong winds have been sug- 
gested also as potential sources of high energy particles 
and gamma rays (e.g. Il4ll5ll6l l29ll30ll . In some of these 
works particle acceleration occurs mainly in the end ter- 
minal shock and in others in shocks embedded in the 
winds. Some of these works predicted detections of iso- 
lated massive stars by EGRET, something that did not 
happen. Efficient electron and proton acceleration in the 
inner wind are unlikely, because of the strong losses to 
which both species are exposed during the processes. In 
addition, even if some gamma rays are produced they 
will be absorbed by the stellar photosphere producing 
electron-positron pairs. 

Models based on colliding winds in massive binary 
systems, on the contrary, seem to be far more promis- 
ing JTl |. The clear detection of non-thermal radio emis- 
sion from the colliding wind region in some systems like 
WR140, WR146, WR 147, and Cyg OB2 #5 indicate that 
electrons are been accelerated up to relativistic energies 
in these binaries. These electrons are located relatively 
close to the secondary star, so they should have strong in- 
verse Compton losses which might result in a detectable 
gamma-ray signal 11321 [33 1. Protons can be in principle 
also accelerate in the colliding wind shock and they can 
interact either with the wind material [34| or they can 
diffuse to interact with nearby molecular clouds ll35l . De- 
tailed models have been recently developed for systems 
like WR140 |[34l[36l . on the basis of the recent radio in- 
formation obtained for this source (37| (see Fig. [2|, and 
for WR20a |38|. This latter system, however, is not re- 
solved at radio wavelengths and since it is very compact, 
strong high-energy photon absorption should take place. 

It is widely expected that GLAST-Fermi satellite and 
even perhaps AGILE could detect some colliding wind 
binaries. Depending on the orbital parameters, these 
sources might appear as highly variable at MeV-GeV en- 
ergies. It is however dubious whether they can accelerate 
particles up to TeV energies. 

GAMMA RAYS FROM MASSIVE 
PROTOSTARS 

In very young SFR, the gravitational contraction of dense 
molecular cores can give rise to a very rich phenomenol- 
ogy. Molecular clouds are clumpy and highly inhomo- 
geneous. In over-dense regions massive protostars are 
formed and ultra-compact regions of ionized gas sur- 
rounds them. The size of these regions is very small, 
< 0.05 pc, implying short lifetimes of the order of 10 3 ~ 4 
yr. The cloud around the protostar is heated and radiates 




FIGURE 2. Evolution of the colliding wind region at radio 
wavelength in WR140 (37). 



in the IR with extremely high luminosities that can reach 
10 5 Lq in some cases. The prestellar core is likely to have 
angular momentum resulting in the formation of proto- 
stellar disks. Evidence for the existence of such disks 
in been found through IR, optical and line observations 
(see, eg., Fig . [3]> . Magnetic fields, which are thought to be 
strong in the inner cores of the molecular clouds (39], are 
dragged in and twisted by the infalling matter. The result 
is a magnetic tower that can launch and collimate a ther- 
mal jet, with terminal velocities of up to a few thousand 
km s _1 . Such outflows are observed in numerous clouds 
(see Fig. HI. The accretion rate at this phase can be as 
high as several times 10~ 3 M Q yr -1 . The whole forma- 
tion process, which is illustrated in Fig.[5]can last ~ 10 5 
yr. 

As time goes by, the star starts to burn hydrogen and 
the subsequent wind blows up the circumstellar disk and 
the jet is quenched. Maser emission is a good tracer 
of the different aspects of the formation process. Three 
molecules are of particular importance: H2O, OH, and 
CH3OH. This maser emission not only acts as a signspot 
of active SFRs, but also provides unique information 
on the physical conditions in the inner regions of the 
molecular core. 

Water maser emission is usually associated with the 
presence of molecular outflows and disks. The lines 
present a wide spread in velocities, that indicating the 
motions of the molecular material in the region. Hy- 
droxyl masers are indicators of the presence of ultra- 
compact HII regions. They are always associated with 
the densest parts of the cloud, being originated in ex- 
tremely small regions of less than 10 16 cm. Finally, 
methanol lines arise in shocked gas that is located at 



the interface between the mass outflow and the external 
medium. Altogether, these maser lines can shed light on 
many aspects of the massive star forming process [ 14|. 

The jets of massive young stellar objects propagate 
through the molecular cloud and can break out from its 
boundaries or can end inside it, remaining hidden at opti- 
cal wavelength. Radio observations, however, reveal the 
presence of jets as thermal, elongated features (e.g. |40|). 
At the end points of the jets strong shocks are expected 
to form. These shocks can be sites of diffusive shock ac- 
celeration, as indicated by the detection of non-thermal 
emission of synchrotron origin in a few sources, like 
IRAS 16547-4247 gHEEl, Serpens (43), W3(OH) flH 
and HH 80-81 BOH . The first case is particularly interest- 
ing, since it seems to be a massive protostar completely 
embedded in a dense molecular cloud with an extreme IR 
luminosity of L ~ 6.2 x 1O 4 L 2.4 x 10 38 erg s" 1 . The 
source is located at 2.9 kpc 0411 . The total mass of the 
cloud is M c \ = 9 x 10 2 M & . Molecular line observations 
indicate that the size of the cloud is ~ 0.38 pc in diam- 
eter (« 1.1 x 10 18 cm). Assuming a spherical geometry, 
the averaged particle (atoms of H) density of the cloud is 
n c i « 5.2 x 10 5 cm~ 3 . The energy density of IR photons 
in the cloud, assuming an homogeneous distribution, is 
w p h w 1.8 x 10~ 9 erg cm~ 3 . 

Radio observations made with the Australia Telescope 
Compact Array and the deeper observations with the 
Very Large Array show the existence of a triple ra- 
dio source inside the molecular cloud (see Fig. [6]). The 
three components of the radio source are aligned in the 
northwest-southeast direction, with the outer lobes sep- 
arated from the core by a projected distance of 0.14 pc. 
The central source is elongated and has a spectral index 
of 0.33 ± 0.05 (S v V a ), consistent with free-free emis- 
sion from a collimated jet [4-2 1 . The integrated emission 
from the northern lobe has a spectral index of —0.32 ± 
0.29, of dubious thermal/non-thermal nature. The south- 
ern lobe, instead, is clearly non-thermal, with an index 
a = —0.59 ±0.15. The inferred linear size for this lobe 
is sa 10 16 cm. The whole situation is depicted in Fig. [7] 

The presence of relativistic particles in such a com- 
plex medium suggest that high-energy emission is likely 
to occur due to relativistic Bremsstrahlung and inverse 
Compton interaction with IR photons in th case of elec- 
trons, and through inelastic pp collisions in the case of 
nuclei. The magnetic field in the region of non-thermal 
radio emission can be estimated assuming equipartition 
and the obtained values are in agreement with direct mea- 
surements by Zeeman effect [39 1. Typical fields are of the 
order of 10 3 G. The acceleration efficiency through diffu- 
sive shock acceleration at the terminal jet region in the jet 
material is small (~ 10~ 5 ) but enough to yield maximum 
energies for primary electrons and protons of the order of 
3 TeV and 60 TeV, respectively, in the case of a proton to 
electron density ratio of 100, similar to what is observed 
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FIGURE 3. A protostellar disk and the associated jet as observed by the Hubble Space Telescope. Courtesy: NASA. 



in the galactic cosmic ray spectrum. The spectral energy 
distribution produced by these particles through the dif- 
ferent cooling radiative mechanisms is shown in Fig. [8] 
In this spectrum not only the contribution from primary 
particles is shown, but also the emission from secondary 
pairs resulting from charged pion decay chains B31 . 

From Fig.[8]it can be seen that integrated luminosities 
up to 10 33 erg s _1 can result from massive protostars with 
strong outflows in a SFR. Such emission can be detected 
with GLAST or with the next generation of Cherenkov 
imaging telescope arrays, opening a new window to the 
study of star formation processes and the environment of 
massive your stars. The combined effect of many proto- 
stars in the core of a very young OB association might be 
even responsible for some unidentified EGRET sources 
found in the plane in positional coincidence with such as- 
sociations [3 1. The higher angular resolution of GLAST 
can differentiate the different individual contributions. 

CONCLUSIONS 

Star-forming regions can generate gamma-rays through 
different processes during their evolution from dark 



clouds to open clusters. In the early stages of star for- 
mation, the contribution from individual protostars can 
reach integrated luminosities above 100 MeV of around 
10 33 erg s _1 . The composed emission of several of such 
YSO can illuminate the massive dark cloud in gamma- 
rays producing a detectable source, even for instruments 
like EGRET and AGILE. Higher resolution telescopes 
like LAT of the GLAST-Fermi satellite might resolve 
the individual sources, then opening a new observational 
window to the study of massive star formation. Strong 
IR luminosities, non-thermal radio emission and maser 
activity are the best tracers for such young protostellar 
objects. 

Once the SFR has aged and early-type massive stars 
are formed, colliding wind binaries can be significant 
gamma-ray sources for GLAST-Fermi satellite or even 
for AGILE. Luminosities up to 10 34 erg s _1 at E > 100 
MeV can be expected for WR-OB systems like WR 140. 
The best indication of efficient particle acceleration in 
such objects is the clear detection of no-thermal radio 
emission from the colliding wind region. Protons might 
escape the system illuminating nearby molecular clouds. 

When the SFR has evolved into an open cluster, the 
collective effect of stellar winds of the strongest O stars 
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FIGURE 4. 
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Star forming region in Carina, with several molecular cores and associated outflows. Image from the HST. Courtesy: 
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FIGURE 5. Stages in the formation of a massive star through 
accretion in a dense molecular core. 



and supernova remnants can produce large-scale effects, 
including particle acceleration up to TeV energies. Such 



heterogeneous sources can be detected in gamma rays 
with Cherenkov telescope arrays as extended sources 
positionally coincident with OB associations and open 
clusters. 

Regarding isolated massive stars, it seems very dubi- 
ous to the author whether they can produce a significant 
amount of high-energy emission. Particle acceleration in 
the inner stellar winds could hardly be efficient, and in 
any case losses, both radiative and adiabatic, make very 
difficult for the particles to reach high-energies. 

Summing up, SFR are promising sites for the produc- 
tion of gamma-rays, by different mechanisms, along their 
whole evolution. In the next few month and years, star 
formation studies might incorporate gamma-ray astron- 
omy to probe the most extreme consequences of the for- 
mation and evolution of massive stars. 
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